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Abstract 

The effect of a uniform dilatation of space on stochastically driven nonlinear field theories is 
examined. This theoretical question serves as a model problem for examining the properties of 
nonlinear field theories embedded in expanding Euclidean Friedmann-Lemaitre-Robertson- Walker 
metrics in the context of cosmology, as well as different systems in the disciplines of statistical 
mechanics and condensed matter physics. Field theories are characterized by the speed at which 
they propagate correlations within themselves. We show that for linear field theories correlations 
stop propagating if and only if the speed at which the space dilates is higher than the speed at 
which correlations propagate. The situation is in general different for nonlinear field theories. In 
this case correlations might stop propagating even if the velocity at which space dilates is lower 
than the velocity at which correlations propagate. In particular, these results imply that it is not 
possible to characterize the dynamics of a nonlinear field theory in an inflationary scenario a priori. 

PACS numbers: 03.50.-z, 64.60.Ht, 89.75.Da, 98.80.Cq 
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s we are going to consider 

HQ. 

In this context one 



Inflation is a cosmological concept that refers to the exponentially fast growth of the 
universe at its early stages. This way, inflation would be responsible for the current homoge- 
neous and isotropic appearance of the universe. In this work we will use the term in a broader 
sense. We will consider stochastic field theories whose correlation length can be explicitly 
computed. When these theories are embedded in an universe which is uniformly expanding in 
time, what mathematically reduces to considering a Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metric, an effective loss of correlation takes place whenever the expansion is fast 
enough. Intuitively, determining the threshold of the speed of expansion that causes this 
loss of correlation would mean to find the condition that ensures the macroscopic appear- 
ance of the universe is homogeneous and isotropic. Although in reality this condition is 
just a necessary but not sufficient one, it is still one of the key ingredients in the search for 
determining the macroscopic appearance of one such universe. This problem has also a very 
natural statistical mechanical motivation. First of all, the mode 
are paradigmatic in the theory of dynamic critical phenomena 
would like to determine whether or not fluctuations are able to break the homogeneity of a 
space which undergoes a uniform dilatation, given that this mechanism acts as a homoge- 
nization on the large scale. On one hand, this question is clearly reminiscent of symmetry 
breaking problems in the context of nonequilibrium phase transitions. On the other hand, 
biological systems have a connection with it too, mainly in the context of pattern formation 
and generation of form during growth si, 0]. Indeed, the equations this work is focused 
on can be considered as hydrodynamic descriptions of the model Eden introduced as an 
idealized description of a growing cell colony p, [6]. 

In this work we are concerned with stochastic fields <f)(x, t) for which the spatial coordinate 
is vectorial x e M. d and the temporal coordinate is scalar t e M + ; in fact we will consider 
an origin of time i > so t is necessarily equal or greater than t . This field will obey the 
generic equation of motion 

d t <f> = f{L 1( f>,L 2( f>,---) + ax,t), (1) 

where the Li's, i = 1, 2, • • • , are linear differential or integro-differential operators acting on 
the field, / is an in general nonlinear function of its arguments and £ denotes a space-time 
noise to be specified in the following. 

Field theories can be characterized by sets of exponents. One of them is the dynamic 
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exponent z which characterizes the velocity at which correlations propagate. If t is the 
correlation length of one such theory then t{t) ~ t x l z \ therefore z > so correlations 
propagate as time evolves. 

We will consider inflation as a simple transformation of the spatial coordinates x — > b(t)x, 
where b(t) is a function of the temporal variable only such that b(t) > 1 for t > t an d 
b(to) = 1. So this transformation is a strict dilatation. 

Field theories can be characterized by different types of correlation functions, such as 
the two-point function and the field difference correlation function. For the models under 
consideration both will undergo dynamic scale invariance; so explicitly they read, the first 
one 

G 2 (x, x>- tf := («/>(*, t)cj>{x>, t)) = \x- xT? (^7^) , (2) 
and the second one 

G d (x, x'; tf := ( \<f>(x, t) - cf>(x', t)\ 2 ) = \x - x'\ 2a Q {^f^J , (3) 

where role of the dynamic exponent is evident and the a exponent describes the variation 
of the field on a determined length scale Q]; J 7 , Q are the scaling functions. 

We say that a linear superposition principle holds whenever, upon applying the inflation 
transformation {x, x'} — > b(t){x,x'}, the correlation functions read 

g\ 24} = b(tf«\x - x'nr,g} ( m *~ x ' 1 ) , (4) 

for the same exponents a and z. It is clear where this expression comes from: in this case 
the internal dynamics of the field and the inflationary effect are simply superposed. 

We start our discussion with the following family of linear equations of motion for the 
field 0: 

d t <f> = -z/|V| c + /i(t)+£(x,t), (5) 
where the noise is assumed to be Gaussian and white, with zero mean and correlation 

(ax,m(x',t')) = D5(x-x f )5(t-t'), (6) 

and the operator IV^ (we will always consider ( > 0) is to be interpreted in the Fourier 
transform sense (|V|^0)" = \k\^<t>. This operator accounts for the anomalous diffusion of 
the field, and its effect on this type of theories has already been considered, even in the 
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nonlinear case [8|. First we note that the appearance of the function h(t) is trivial as it can 
be adsorbed by means of the transformation <\> — > <\> + J h(t)dt and we will do so in the 
following. This linear field theory can be exactly integrated and shown to obey the above 
mentioned scalings with z = C (note that for ( integer and even the operator becomes, up 
to its sign, the Laplacian or a power of it). Actually, the presence of noise in this equation 
is trivial in the sense that the exponent z does not change if we set D = (trivial in this 
respect, not necessarily in others). 

Now we will apply the inflation transformation to Eq. (jSJ), and for the sake of concreteness 
we set b(t) = (t/t ) 7 , for 7 > 0, and we will refer to this parameter as the growth index. The 
natural question that arises is: is there linear superposition for this case? Even in this case 
in which the equation is linear the answer is only partially positive (and so in general it is 
negative). Linear superposition holds for both correlation functions, without amendments, 
only if 7 < 1/0 For 7 > 1/C the two-point function does not adopt the form described 



by Eq. (j4j) [4j, Ildl-ll2|. Indeed, this value of 7 plays a special role. For 7 < 1/C one can 
read from Eq. (jlj) that correlations still propagate as time evolves. The contrary would 
happen if we reversed the inequality. So the linear superposition principle has a transparent 
physical meaning: propagation of correlations and dilatation of space are simply superposed. 
Although linear superposition does not take place for the two point function and large enough 
7, we still have a weaker yet intuitive result: propagation of correlations stops (or at most 



becomes logarithmically slow) in the limit 7 — > (1/C) generically in the general case 



we should say that propagation of correlations stops whenever 7 > as some particular 
cases admit a different marginal propagation of correlations for 7 = 1/C {12I . [13I ] . 

Explicitly, after applying the inflation transformation, the equation of motion reads 

M = ~ (7) ' 7 + (7) d1 ' 2 t). (7) 

Following Eq. (jlj) and for 7 < 1/C one may define an effective dynamic exponent z e g = 
C/(l — 7C) [Uj- Thus, in particular, in the limit 7 — > + one recovers the classical case 
z e ff — > C, and when 7 — > (1/C) - then z e s — >■ 00. So we may talk about the decorrelation 
threshold 7^ = 1/C- As we have already mentioned, although linear superposition does not 
completely hold in this case, the decorrelation threshold is the same as if it held. 

The main conclusion of this analysis is that for a rather general family of linear equations 
the decorrelation threshold is the intuitive one. We will show that for nonlinear equations 



things are in general different. To this end one needs to introduce a nonlinear field theory 
whose dynamic exponent is nontrivial. One such theory is given by the Kardar-Parisi-Zhang 
(KPZ) equation |l| 

^ = z/V 2 + ^(V0) 2 + e(x,t). (8) 
Together with the interest of this equation in the fields of condensed matter and statistical 



physics 



151 ]. one finds its relevance in cosmology and related topics 



16-241. Two of the rea- 



sons underlying this universal character are the connection of Eq. ([8]) with Burgers equation 
through the definition of the velocity field v := and to the imaginary time Schrodinger 
equation with a random potential by means of the change of variables ip := exp[A0/(2z/)]. 
If we set D = in this equation we find that z = 2, as can be read from its exact solu- 
tion in this case [yj], or even by simple dimensional analysis. However, contrary to what 
happened to the linear equations, once the noise is switched on, the dynamic exponent 
becomes a function of the spatial dimension, z = z(d). In particular, it is known that 
z(l) = 3/2 and z(2) ~ 1.7 < 2. We note that while the one- dimensional result is exact, 
the two-dimensional one is usually obtained with numerical methods. This is so because 
the calculation of this value has escaped all sorts of analytical approaches, with the notable 
exception of an improved-renormalization-group-type method known as the Self-Consistent 



Expansion (SCE) 25j, |26j. Due to the remarkable success of this scheme in finding the 



scaling behavior of this as well as different models 27H29| we will rely on its results in the 
following. 

Our aim is calculating the decorrelation threshold for the KPZ equation. Following the 
linear theory one could nai'fly expect 7^ = 1/z, and in particular — 2/3 in d — 1. In fact, 
the one dimensional result is presumably correct. Simulations of a discrete model in the 
KPZ universality class have corroborated this result {3^]. Things are, on the other hand, 
different in higher dimensions as will shown in the forthcoming analysis. 

We now apply to Eq. ^ the dilatation transformation x — > (t/t o y x, where as always 
7 > 0. The resulting equation reads 

+ \ 2 7 x / . \ 2 7 , , x d 7 /2 

to \ 1 1 / \ /v7 A\2 1 / r 



dt4> = »[j) ^+2^7 J t(x,t). (9) 

So we will study this equation which describes the KPZ dynamics in an environment which is 
undergoing spatial dilatation as time evolves. This is a Langevin equation whose associated 
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Fokker-Planck equation reads 



ihV = ( ^ 



2 7 



dx— 



^V 2 0+^(V0) 2 



dx-—V, (10) 



where the solution V is the functional probability distribution. This equation can be trans- 
formed to 



dV 



I 



dx 



5(f) 



z/V 2 0+^(V0) 2 



VM D ft V d - 2h f , 5 2 



(ii) 



Now we change the temporal variable 

27 

dr — I — 1 (it 



27 



> r - j^^t t J- 

We start assuming 7 < 1/2 so that r is growing for growing t and 



(12) 



t 27 

r ps — - — t 1 27 when t ->■ 00. 
1 - 27 



(13) 



After performing the change of variables and going back to the Langevin equation we find 

(14) 

* \ L o / 

or alternatively 



(2-d) 7 

A „ / 1 — 2^y\ 2(1-2^) (2 -d) 7 

,9 T = z/V 2 + ^(V0) 2 + r2TT^)^(x,r) 
2 V *o / 



d T = z/ V 2 + ^(V0) 2 + r ? (a;,r), 



where the noise correlations are (t](x,t)) = and 



(2-d)7 



(77(2;, t)t7(x', r')) = D t ^ 5(x - x')S(t - r'), 



(15) 



(16) 



for a suitable constant D. This noise rends in general a more difficult analytical treat- 
ment due to the explicit r— dependence of its amplitude. However, the situation becomes 
considerably simpler in d = 2. In this case we have 



9 T = z/V 2 0+^(V0) 2 + e(a:,r) ) 



(17) 



which is the KPZ equation for time r, so the dynamic exponent for this time variable is 
z' = ZKPz(d = 2). Thus the effective dynamic exponent for actual time t is 

ZKPz(d = 2) 



ZeS 



1 - 2 7 



(18) 
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Consequently decorrelation appears in the limit 7 — >■ (l/2)~ =>■ 7^ = 1/2, instead of 
7d = l/zKPz(d — 2) « 0.6 > 1/2. So the decorrelation threshold is anticipated, and 
this counterintuitive result implies that a simple superposition principle does not hold in 
this case. 

Although this result proves the decorrelation threshold by itself, it is easy to compute the 
exact behavior at the value 7 = 1/2. In this case one can correspondingly modify change 
of variables f fT2|) to find r = £0 hi(£/£o)- So for the critical value of 7 correlations propagate 
logarithmically slow, and we find the effective value z e ff = 00. For 7 > 1/2 change of 
variables (IT2"|) is still valid. One again finds that the solution of Eq. OH]) becomes the solution 
of the classical KPZ equation in time r. The particularity of this situation is that, as time 
£ progresses, time r evolves from r = when £ — £q to the finite value r = £0/(27 — 1) in 
the limit £ — > 00. So the resulting profile of the solution to Eq. (Q becomes the profile of 
the solution to Eq. ((§]) quenched at time £0/(27 — 1) asymptotically in time. This result, 
among other things, implies that correlations propagate for short times but tend to stop 



asymptotically in time. This means that the di 



atation of space acts as an effective drift 



that is able to transport correlations by itself jll ]. 

This effect is not purely two-dimensional. We now move to higher dimensions and consider 
again the KPZ equation but with a different stochastic forcing 

^ = ^V 2 0+^(V0) 2 + X (M), (19) 

where the noise is assumed to be Gaussian, white in time and spatially correlated. We 
consider it has zero mean and explicitly its correlation reads 

( X (x,t)x(x',t'))=D\x-xfP- d 5(t-t'), (20) 

where p > specifies the degree of spatial correlation (p = sends us back to the spatially 
uncorrelated noise). This equation can again be mapped onto a Fokker-Planck description. 
And again, the same transformation r = £ 2 1 7 (£ 1-27 — £q~ 27 )/(1 — 27), yields in the limit £ — > 00 
the equation 

(2+2p-d) 7 

A, x0 fl — 27 \ 2 ( 1 - 2 ^) (2+2p-d) 7 



<9 r = z/V 2 + -(V0) 2 + [—^) T^^r x (x,T). (21) 

In this case this model becomes exactly Eq. ( TT9|) in time r for d = 2 + 2p. So choosing 
appropriate integer and semi-integer values of p one recovers the KPZ equation in r time 



for any desired spatial dimension d > 2. Model (fT9|) was analyzed with the SCE and for 
the dimension under examination d = 2 + 2p classical KPZ behavior, as if p = 0, was 



found 



311 ] . It is not clear whether or not there exist an upper critical dimension for KPZ 



(a dimension above which the large-scale effective behavior of the equation wou 
that of its linear counterpart) and what would be its value in the first case 32 



d reduce to 



331 ] . Recent 



numerical results suggest that, if it exist, one necessarily has d c > 4 [34j. In any case, it 
is clear that for any d > 2 under the upper critical dimension of KPZ the corresponding 
dynamic exponent of model (TT9|) is z < 2, while the decorrelation threshold is as before 
anticipated and results — 1/2. This proves that the nontrivial coupling of the inflationary 
transformation and the nonlinear field theory extends from two to higher dimensions, at least 



in the range 2 < d < 4 according to 



341 ]. and possibly to higher dimensions. This also shows 



that apparently the one-dimensional situation is left alone as the only one in which the 
decorrelation threshold is the intuitive one. And thus, this fact is yet another proof of the 
fundamentally different character of the KPZ equation in and above one dimension, posing, 
in the latter case, a problem much more involved and changeling. This difference should 
be even more pronounced in the neighboring field of radial growth, since posing the KPZ 



12 



problem in this context implies nontrivial topological effects only if d > 2 

As already mentioned in the introduction, the condition 7 > 1/z implies correlations 
stop propagating in the linear case, but it is just a necessary and not sufficient condition to 



achieve the spatial homogeneity of the fie 
if 7 > max{l/z, 1/d} in the linear case 



?ld. Homogeneity is only achieved in the large scale 
4j, what shows that the spatial dimensionality of 
the system has an important role in this question. Although it would be tempting to simply 
extend the previous results concerning the linear models to the general case, we are going 
to show that such an extension is not correct. However, we will be able to find an analogous 
condition by means of the introduction of new critical exponents a and z; furthermore, these 
new exponents will allow expressing the condition for decorrelation as 7 > 1/z. To this end 
we will make explicit use of correlations (jSJ) and (]3]) for the KPZ case. If we write these 
correlations in the form suggested by Eq. (jl]), so that the dependence on the dilatation of 
space becomes explicit, we find the expression 

't 7 |x — x'\ 



r 2 



t 2a 7 | 



J\2a 



6} 



t l/z 



(22) 



where the new exponent values are a = (1 — 27)«kpz/ (1 — 27 + z^pzl) and z = zkpz/ (1 — 
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27 + zkpzi) whenever 7 < 1/2, where «kpz and zkpz are the corresponding exponents of 
KPZ at the lower critical dimension d = 2 + 2p (including the case p = 0). For 7 > 1/2 
the exponents become at = and 5 = 2, with marginal logarithmic corrections for 7 = 1/2. 
Also, if one uses the form of both correlations G = t 2/3 H(|x — x'\/t l l z ) one may extract 
in our case the value of the new growth exponent $ = (1 — 27)/3kpz for 7 < 1/2 and 
$ = for 7 > 1/2 (with again a marginal logarithmic correction for 7 = 1/2) and find 
the relation — 6l/z holds in this case too. Additionally one sees that using the relation 
«kpz + zkpz = 2 one finds for the new exponents at + z = 2. The different relation found for 



the KPZ equation with a time dependent coefficient of the nonlinearity 35| can be written 
in the present terms as «kpz + z cS = 2 + 2jz e g. From these results one reads that the 
dilatation of space systematically neglects the non-perturbational behavior of KPZ. A new 
quantity within reach is the center of mass fluctuations, which can be characterized by a new 
exponent J 4>(x, t)dx ~ t M (this time the integral should be carried out on a bounded domain 
for obvious reasons). In the present case this exponent reads fi — \J3+ (d — 2p)/(2z)](l — 27) 
for 7 < 1/2 and p = for 7 > 1/2, with as previously a logarithmic correction for 7 = 1/2. 
In particular, the motion of the center of mass becomes bounded in time for 7 > 1/2. Note 
that this quantity is closely related to the properties characterizing weak convergence of the 
profile of (p to the homogeneous spatial state |4| . In general the exponents which characterize 
weak convergence to the homogeneous profile are different from the exponents appearing in 
the field difference correlation function ([3]) [4J. However, in the present cases, both sets of 
exponents are exactly the same. This should be considered by no means a general feature 
of the KPZ equation: it is a direct consequence of fact that we are always considering this 
equation at its lower critical dimension. Note also that our results imply the flatness of the 
field (in the sense that both field difference and two-point correlation functions are uniformly 
bounded in both space and time) is achieved for 7 > 1/2. This is in disagreement with the 
linear requirement 7 > max{l/z, 1/d}. However, it would be in perfect agreement with the 
modified requirement 7 > max{l/z,l/d}. In the same way, the threshold for the loss of 
correlation could be expressed by the inequality 7 > 1/5. Both inequalities express the fact 
that the linear conditions for decorrelation and homogeneity of the field can be extended to 
the nonlinear case provided we introduce 7— dependent exponents. This is another way of 
expressing that the coupling of the inflation transformation and the nonlinear dynamics of 
the field is nontrivial. 
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We also note that, although we have proven the threshold 7 = 1/2 for decorrelation 
and homogeneity of the field for the lower critical dimension of KPZ, we expect it will stay 
the same for dimensions above this one. This conjecture comes from the fact that change 
of variables f[T2"j) sends the equations under consideration to KPZ equations with weaker 
noises (noises whose amplitude depend on a negative power of time) in the case of a higher 
dimensionality. The question of super-roughness of the field, i. e. finding the values of 7 
for which the fluctuations of the field grow faster than the dilatation of space, is a simple 
corollary of our results. This would happen whenever a > 1, what is impossible for any 
7 > 0. 

In summary, we have studied the effect of a uniform dilatation of space, what we have refer 
to as the inflation transformation, on the dynamics of nonlinear fields theories. In particular 
we have focused on the nonlinear KPZ equation with different stochastic forcing terms, 
because this field theory is known to display nontrivial effects regarding the velocity at which 
correlations propagate. We have argued that in one dimension numerical results suggest 
that the loss of correlation starts when the velocity at which the space grows overtakes the 
velocity at which correlations propagate in the absence of inflation. However, in two and 
higher dimensions the threshold for the appearance of decorrelation becomes anticipated, 
and so loss of correlation starts at a velocity of the inflation transformation slower than the 
speed at which correlations propagate. This fact is a consequence of the nontrivial behavior 
of the KPZ equation at the lower critical dimension. It shows that the interplay of inflation 
and nonlinearity is far from trivial and, in particular, that it is not possible to infer the 
effect of a dilatation of space on a nonlinear field theory a priori. 

There are several interesting connections among models in condensed matter physics and 
cosmology. In this work we have discussed one such model given by the KPZ equation, 

and Bose- 



which lies in the mentioned interface as well as Ginzburg-Landau theories | 
Einstein condensation 37]. Some questions naturally emerge from the present study. One 
is determining under which conditions loss of correlation in an anisotropically expanding 
universe is achieved. Mathematically, accounting for anisotropic expansions implies the 



substitution of the FLRW metric by a Bianchi I metric 



38] . Another problem is the analysis 



of related nonlinear models with a source of quantum fluctuations instead of the classical 



ones. In this framework the question of under which conditions disentanglement occurs 
seems to be connected with the present discussion. 



39] 
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